Simple quantum error detection and correction for superconducting qubits 
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We analyze simple quantum error detection and quantum error correction protocols relevant to 
current experiments with superconducting qubits. We show that for qubits with energy relaxation 
the repetitive iV-qubit codes cannot be used for quantum error correction, but can be used for 
quantum error detection. In the latter case it is sufficient to use only two qubits for the encod- 
ing. In the analysis we demonstrate a useful technique of unraveling the qubit energy relaxation 
into "relaxation" and "no relaxation" scenarios. Also, we propose and numerically analyze several 
two-qubit algorithms for quantum error detection/correction, which can be readily realized at the 
present-day level of the phase qubit technology. 
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I. INTRODUCTION 

Quantum error correction [l[ (QEC) is an unavoidable 
procedure in a practical quantum computer [2, 3]. The 
standard QEC 0, includes encoding a logical qubit 
in several physical qubits, measuring the error syndrome 
using ancillary qubits, and then applying a correction op- 
eration, which depends on the measurement result. (A 
promising variation of this idea are the so-called surface 
codes [J]-) Unfortunately, QEC is very difficult experi- 
mentally [5l4l5l| , and some simplifications are often used. 
Let us mention three of them, all of which have been in- 
troduced in Ref. First, instead of using additional 
qubits for the error syndrome, in a "compact" scheme 
the same physical qubits can be used for the encoding 
and error syndrome measurement; this is done by de- 
coding the encoded state after a possible error occurs. 
Second, since a single-shot measurement of a qubit state 
is often difficult, the standard QEC can be replaced by 
measurement-free QEC, in which the measurement and 
correction are substituted by a quantum controlled op- 
eration (e.g. the Toffoli gate). Third, a favorable type of 
error (against which the code protects) is often simulated 
by applying a certain unitary rotation, with the rotation 
angle corresponding to the error strength. 

Measurement-free QEC experiments in nuclear mag- 
netic resonance (NMR) systems [^-Q have been per- 
formed for over a decade, but only with ensembles of 
quantum systems 16]. Using trapped ions, a three- 
qubit QEC experiment with actual measurement was re- 
alized [9]. and recently a measurement-free QEC pro- 
cedure with several error correction cycles was demon- 
strated In linear optics systems, the QEC exper- 
iments include t wo-q ubit protection against "acciden- 
tal" measurement [ll| , continuous- variable adaptation of 
the 9-qubit Shor's code [12j . continuous- variable erasure- 
correcting code [l3j], and eight-photon topological error 
correction flij , A three-qubit measurement-free QEC 
protocol has been recently demonstrated with supercon- 
ducting "transmon" qubits [iEj . 

With the rapid progress in experiments with supercon- 



ducting qubits [13, [l8[ , QEC with actual measurements 
is becoming feasible in these systems in the reasonably 
near future. The subject of this paper is the analysis of 
several simple quantum error correction/detection proto- 
cols relevant to future experiments with superconducting 
qubits, mainly superconducting phase qubits 1191. ( Some 
results of this paper have been reported earlier |2fj|.) 

In the past, pure dephasing was by far the dominant 
source of decoherence in superconducting qubits, and 
QEC protecting against pure dephasing would be most 
important. An example of such a procedure was consid- 
ered theoretically in Ref. [2l|. The idea was to use the 
standard 3-qubit repetitive code, which protects from bit 
flips (i.e. X-rotations). By using additional Hadamard 
gates for each physical qubit, the X-rotations are con- 
verted into Z-rotations, and therefore the same code can 
be used to protect against pure dephasing. 

In recent years, pure dephasing in superconducting 
qubits was si gnif icantly reduced by various technological 
advances [i~7lTl8j , and now energy relaxation is becoming 
most important. In particular, when quantum informa- 
tion is stored in a superconducting resonator [22|, l23j . 
pure dephasing is negligible in comparison with energy 
relaxation. This is why in the first part of this paper (Sec. 
II) we focus on the operation of repetitive iV-qubit quan- 
tum codes in the presence of energy relaxation. Repeti- 
tive codes are chosen because of their relative simplicity 
in the encoding and decoding (unfortunately, the stan- 
dard 5-qubit or 7-qubit stabilizer codes 0, [H, [24| are not 
feasible for superconducting qubits in the near future). 
To reduce the number of qubits in a procedure we use 
the standard compact scheme [f| [2lJ, in which the an- 
cilla qubits used for encoding are also used for the er- 
ror syndrome measurement. We assume that the energy 
relaxation happens at zero temperature, which is essen- 
tially the case for superconducting phase qubits, since 
the typical qubit frequency is ~6 GHz, and therefore the 
energy Tiu) ~ 0.3 K is much larger than the experimental 
temperature of ~50 mK. 

Even though energy relaxation may look similar to a 
bit-flip, it actually can be thought of as a combination 
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of two quantum errors: bit- flip and bit-phasc-flip (which 
correspond to A-rotation and Y- rotation). This is the 
reason why, as we show later explicitly, repetitive codes 
do not work for QEC against energy relaxation. How- 
ever, these codes can be efficiently used for quantum er- 
ror detection (QED). In QED we detect that an error 
happened but cannot restore the undamaged quantum 
state (in particular, the QED idea was implemented in 
Ref. [25| for phase errors in liquid-state NMR and has 
been recently investigated in Ref. [l3| for detecting pho- 
ton erasures). Even though QED is of much more limited 
use than QEC, it is still an interesting procedure, and ex- 
perimentally can be considered as a first step towards full 
QEC. We show that for QED against energy relaxation 
it is sufficient to use 2-qubit encoding and that there is 
not much benefit to use more qubits, unless a somewhat 
more sophisticated procedure is used. 

Our analysis in Sec. II is based on unraveling the qubit 
energy relaxation into the "relaxation" and "no relax- 
ation" scenarios. This unraveling is quite different (and 
more difficult) than, for example, unraveling of pure de- 
phasing into the "phase flip" and "no phase flip" sce- 
narios. The main reason for the difference is that the 
unraveled states for the energy relaxation are related to 
the initial state in a non-unitary way. 

In Sec. Ill we focus on simple two-qubit QEC/QED 
protocols, somewhat similar to those in Ref. [25|, which 
can be readily implemented using the present-day tech- 
nology of phase qubits (H). Realistic experimental pa- 
rameters are used in the numerical simulation of these 
protocols. In the first protocol and its variations, we as- 
sume that as in most of the previous experiments [Ml 
the errors are intentionally induced by particular oper- 
ations pretended to be unknown. The algorithms can 
mainly be used for QED; however, when the type of par- 
ticular error is known (which is the case for intentional 
errors in an experiment), the algorithms can also be used 
for QEC. We also analyze numerically the operation of 
a protocol, in which the errors during a storage period 
are due to actual energy relaxation of two qubits (assum- 
ing storage in resonators of a RczQu-architccturc device 
[23I I2H ] with phase qubits). This protocol can only be 
used in the QED mode. The main result of the simula- 
tions is that the analyzed protocols can be realized at the 
present-day level of phase qubit technology. Section IV is 
a conclusion. Some mathematical details of the analysis 
are discussed in the Appendix. 



II. REPETITIVE CODING FOR ENERGY 
RELAXATION 

In this section we analyze the operation of repetitive 
A-qubit encoding in the presence of (Markovian) zero- 
temperature energy relaxation. The procedure is shown 
in Fig. 1. The goal is to preserve an arbitrary initial state 

|Vi„) =a\O)+0\l) (1) 
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FIG. 1: iV-qubit repetitive coding algorithm with one con- 
trol qubit initially containing the quantum information. The 
controlled- V block represents CNOT gates from the main 
qubit to each ancilla qubit individually. represents en- 

ergy relaxation of the ith qubit (i = 1 for the main qubit, 
i > 2 for ancilla qubits). 



of the main (upper) qubit, where |0) is the ground state 
and |1) is the excited state. In this paper we consider only 
preservation of the initial state ( "memory" operation) , so 
in discussing the fidelity of a procedure we always imply 
comparison with the ideal memory operation. 

The encoding in Fig. 1 is performed with N — 1 
controlled-NOT (CNOT) gates, acting on N - 1 ancilla 
qubits, which all start in the state |0). This produces 
the iV-qubit wavefunction a\0 N ) + f3\l N ), where the no- 
tation \x N ) represents the product-state of N qubits, all 
being in the state x. After the encoding, all qubits are 
subjected to decoherence due to zero-tcmpcrature energy 

(i) 

relaxation with relaxation time T± for the ith qubit, 
i = 1, 2, ... N. We will mostly consider the case when the 
decoherence is the same for all qubits, T X W = T x . After 
the decoherence during time t, the logic state is decoded 
by using N—l CNOT gates in the same way as was done 
for the encoding, and all N— 1 ancilla qubits are measured 
in the computational basis. In the absence of decoherence 
(t = 0) the state after decoding is (a\0) + /3|1)) \0 N ~ 1 ), 
so that the initial state of the main qubit is restored and 
the measurement results for all ancillas are 0. The de- 
coherence disturbs the final state, which probabilistically 
changes the measurement results and the corresponding 
final states of the main qubit. 

Even when the measurement result is all N — 1 zeros 
(for which we will use the bold-font notation 0), the state 
of the main qubit is not exactly |f/ , in); however, we will 
see that it is close to \ipin}- A measurement result differ- 
ent from indicates an error. There are three ways to 
handle this situation. First, the measurement result can 
be simply ignored; in this case there is obviously no bene- 
fit from using the encoding/decoding procedure. Second, 
we can reject such cases and keep only realizations with 
the measurement result 0; we will refer to this selective 
procedure as quantum error detection. Third, we can ap- 
ply a quantum operation to the main qubit to make its 
state closer to \ipin}- This operation will depend on the 
measurement result, and the procedure is then quantum 
error correction. 

For simplicity in this section we neglect decoher- 
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ence (and other imperfections) during encoding, decod- 
ing, and measurement; it will be taken into account 
in the next section when we will discuss realistic ex- 
periments with phase qubits. To characterize the effi- 
ciency of a procedure either the quantum process to- 
mography (QPT) fidelity F x or the average state fi- 
delity F av can be used. The QPT fidelity is usually 
defined as [1 F x = Tr(xdcsircdX), where x is the 
process matrix and Xdesired in our case corresponds to 
the ideal quantum memory operation, i.e. no evolution 
of the logic qubit. The average state fidelity is 0, [13] 
-Fav = /Tr(p fin C/o|7/'i n )(-0i n |f7 n t )d|i/'i n ), where U = 1 is 
the desired unitary operator, /Ofi n (|V'in)(^'in|) is the ac- 
tual mapping from the initial state to the final density 
matrix pfi n , and the normalized integral is over all pure 
initial states \ipi n ) using the Haar measure. For a trace- 
preserving operation F av = (F x d+1) / (d+1), where d = 2 
is the dimension of our Hilbert space [27[ ■ This rela- 
tion holds for QEC and/or when the measurement result 
is ignored. However, for the QED procedure there is 
a problem [29] in defining the QPT fidelity F x because 
the procedure is selective; then the quantum operation 
for normalized states is not linear and the corresponding 
(trace-preserving) matrix x cannot be defined rigorously. 
In this case we define F x via the average state fidelity, 

F x = (3F av - l)/2, (2) 

as for a trace-preserving operation. 



A. Single-qubit relaxation 

Before calculating the fidelity of the QEC and QED 
procedures, let us calculate the quantum memory fidelity 
of a single qubit, without any encoding. We also consider 
first this simple case to demonstrate a technique of un- 
raveling the evolution due to energy relaxation, which is 
later used for iV-qubit encoding. 

After time t an initial state \ipi n ) = a|0)+/3|l) becomes 
a density matrix (here the upper row and left column 
correspond to the excited state |1)) 



Pfin 



\2 e -t/ Tl 



a *p e -t/2 Tl 

I 2 + \t3\ 2 (l-e t l T -) 



(3) 



which can be represented using the Kraus operators: 

p&n = A T p in Al + A nPin Al, (4) 
A r 



f 



-4„ 



Vl^P 
1 



(5) 



where p = 1 — e _t ' , p- m = |^m)(^ , in| ) and the Kraus op- 
erators satisfy the completeness relation A\A X + A^A^ = 
1. This representation has an obvious interpretation as 
two scenarios of the evolution. The first term in Eq. 
((3]) corresponds to qubit relaxation into the ground state 
|0) with probability P T = \f3\ 2 p. The second term is the 
no-relaxation scenario, which occurs with the remaining 



probability P n = \a\ 2 + \(3\ 2 (l — p) = 1—P X and transforms 
the qubit into the state 



A n |^i„) a\0) + 0y/T=p\l) 



(6) 



The non-unitary evolution |i/'in) — > is essentially the 
same as for a partial collapse due to a null-result mea- 
surement in the experiment of Ref . [30l | . 

Now let us find the averaged state fidelity F av = 
Tr(pt |^i n )(^in|) using unraveling into the relaxation and 
no-relaxation scenarios. With probability P r the state 
fidelity is F st , r = |(0|i/'in)| 2 = |ct| 2 , and with prob- 
ability P n the state fidelity is F st ^ n = |(V'n|V , m)| 2 = 
(\a\ 2 + y/1 — p\(3\ 2 ) 2 /P n . Therefore for an initial state 
iV'in) the state fidelity is 

F s t = F st>I P x + F stt nP n (7) 
= \a\ 2 \(3\ 2 p+ |«| 4 + (1 -p)|/?| 4 + 2\a\*\/3\ 2 y/T=3, (8) 

and the average fidelity F av = F st can be calculated by 
averaging |cn| 4 , |/3| 4 , and |a| 2 |/3| 2 over the Bloch sphere. 
These averages (including some others for completeness 
and later use) are 



4 - I/PI4 _ 



77 (1 + cos6») 2 sin 9 1 
da = — . 



(9) 



1 J 4 2 6' v ' 



a- 



2' M 



a- 



l«| 4 |/5| 2 = ^, M 4 | 



1 

30' 



(11) 
(12) 
(13) 



1/(A + B\/3\ 2 ) = (l/B) ln(l + B/A), 
WKA + bW) = (1/-B) - (A/B 2 ) ln(l + B/A), (14) 



a 





1 


A A 2 


l/?l 2 ) = 


2B 


~~B 2 ~ + ~EP 


-3 


A 


(A + B) 2 


" 2B ~ 


B 2 


B 3 
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A 


A{A + B 


= 2B + 


B 2 


B 3 



ln{l + B/A), (16) 
-ln(l + i?M),(17) 



where A and B are constants, and we used integration 
over the Bloch-sphere polar angle 6, so that \a\ 2 = (1 + 
cos6»)/2 and \(i\ 2 = (1 -cos0)/2. 

Applying the averages © and (flUl to Eq. (|5J), we ob- 
tain the average state fidelity 



yi -p _ p 

3 6' 



(18) 



Actually, there is an easier way to obtain this result. In- 
stead of averaging _F st over the Bloch sphere, it is suffi- 
cient [13, HH (see also Appendix) to calculate the aver- 
age only over 6 initial states: |0), |1), (|0) ± |1))/V2, and 
(|0)±i|l))/ v / 2. However, in our further analysis this trick 
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does not always help, so we prefer the full integration over 
the Bloch sphere. Using Eq. @ it is easy to convert Eq. 
(flgf into the QPT fidelity: F x = (1 + V^P ~ P/ 2 )/ 2 - 
Note that for small p 



1 



V - F 

3' x 



f 



2" 



P ■ 



Ti 



< 1. 



(19) 



The average state fidelity (fT5|) is averaged over the two 
scenarios. Let us now discuss the average state fidelity in 
each scenario separately, having in mind a gedanken ex- 
periment in which an emitted photon or phonon is always 
captured and recorded, thus allowing us to distinguish 
the two scenarios. If the relaxation has happened, then 
-Fst.r = |o;| 2 and averaging this over the Bloch sphere we 
obtain 

Fav,r =W= 1/2- (20) 

Similarly, for the no-relaxation scenario -F a v,n = 
(H 2 + i/T=p~\f3\ 2 ) 2 /[\a\ 2 + (1 -p)\/3\ 2 }, which can be 
calculated using Eqs. (fT5 |l - ([T7|l : 

1 , v T^p(2-p)-2(l-p) 



(l-p)(2,/T=p-2+p) 



P 



ln(l-p) (21) 



For p <C 1 this gives P a v,n « 1 — p 2 /24, showing a 
slow, quadratic in time decrease of fidelity in the no- 
relaxation scenario in contrast to the linear decrease (Til?)) 
of the fidelity averaged over both scenarios. Therefore 
our gedanken experiment could be used for quantum er- 
ror detection: if no relaxation is recorded, we know that 
the initial state is well-preserved at short times. 

Note that we have averaged the state fidelities F s t, r and 
Pst,n over the Bloch sphere with uniform weight, as in the 
standard definition @, H3 of the averaged state fidelity. 
Another meaningful averaging is using weights propor- 
tional to the probabilities of the corresponding scenarios. 
[This would correspond to an equal number of experi- 
mental runs for each point of a uniform mesh on the 
Bloch sphere, as opposed to an equal number of "suc- 
cessful" (i.e. selected) runs for the previous definition.] 
Thus defined average fidelities are 



= M 2 P r /p r = 1/3, 



(22) 



^av,„ = (M 2 + V^P~\P\ 2 ) 2 /Pn = 2 3 ^ + 3^P^ (23) 

where P r — p/2 and P n = l—p/2 are the averaged prob- 
abilities of the two scenarios. The advantage of this def- 
inition is a natural formula for the non-selected average 
fidelity: 



-Pav — -Pav.r-Pr ~t~ -Pav,nPn 



(24) 



[see Eq. (UJ]. In this paper when discussing selected sce- 
narios (as for QED) we will use both ways to average 



over the Bloch sphere. Note that F a v,n ~ 1 — p 2 /24 for 
p« 1, which is the same as for P a v,n (F aVi „ and P av ,n 
are practically indistinguishable at p < 1/2), indicating 
that the difference between the two definitions is not very 
significant in the cases that are of most interest for this 
paper. 



B. Two-qubit encoding 

Let us use the procedure of Fig. 1 with only one ancilla 
qubit. The encoded state is then a] 00) + /9|H). The 
state evolution due to energy relaxation can be unraveled 
into four scenarios: no relaxation, relaxation in either 
the first (main) or second (ancilla) qubit, and relaxation 
in both qubits. The corresponding wavefunctions and 
probabilities after time t of energy relaxation are 



( ajoo) + /vr~^yr~^|ii) 

V 1 nn 

prob. P nn = \a\ 2 + \(3\ 2 (l- Pl )(l-p 2 ), 
|01), prob. P m = \(3\ 2 Pl (l-p 2 ), 
|10), prob. P m = \(3\ 2 (I ~ Pl )p 2 , 
I |00), prob. P rr = \P\ 2 p lP2 , 

where 



pi = 1 



-tm L 



> P2 = 1 



-t/Tf 



(25) 
(26) 



are the single-qubit probabilities of relaxation from the 
excited state |1). This simple unraveling is possible be- 
cause the energy relaxation occurs only in component 
1 11) of the superposition, and in this component the 
qubits are unentangled. This is why the probabilities of 
the scenarios are the simple products of individual prob- 
abilities. The validity of Eq. (1251) can also be checked by 
considering particular time moments at which the relax- 
ation events happen and integrating over these moments; 
this is a more direct but more cumbersome way. 

After the decoding procedure consisting of one CNOT 
operation, the two-qubit state is a product-state in all 
four scenarios: 



a|0) +/V^"piVT^|1) 



|01), prob. P rn , 
|11), prob. P nr , 
|00), prob. P rr , 



<g> |0), prob. P n 



(27) 



with a definite result of the ancilla qubit measurement 
in each scenario. The state of the main qubit is dif- 
ferent from the initial state \tpin) in all four scenar- 
ios, and the corresponding state fidelities are [|a| 2 + 

I 2 , and lal 2 . 



2 v /r^ v /r^] 2 /p n 



at 



As discussed above, we consider three ways to proceed: 
ignore the measurement result, select only result 0, or 
try to correct the main qubit state. If the measurement 
result is ignored, then all four scenarios are added up and 
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the average fidelity is 



FT = [|a|* + W^/T^y/T 



P2\ 



H 2 P rl 



a 2 P„ 



(28) 



where the averaging is over the Bloch sphere. Using the 
formulas for the probabilities from Eq. (I2"5t and the aver- 
ages |/3p, jap, j/9j* and |a| 2 |^| 2 from Eqs. ©-(03]), 
we obtain 



F lgn = 



V(l-Pl)(l-P2) 



Pi 

6 



(29) 



For small pi. 2 (at short time t) it is F}& n 1— pi/3— P2/6, 
and it is obviously worse than the case without encod- 
ing/decoding of the main qubit - see Eqs. (fTSj) and ([T9]) . 
Note that Eq. (]2T)]) can also be obtained by averaging the 
state fidelity only over the 6 initial states (see Appendix). 

In quantum error detection we consider ancilla mea- 
surement result 1 as an error and select only the cases 
when the measurement gives 0. This selects scenarios 
with either no relaxation or two relaxation events [see 
the first and last lines of Eq. (|27|) ]. The averaged (with 
uniform weight) state fidelity in this case is 



av P + P 

J nn J rr 



(30) 



(the fraction is averaged over the Bloch sphere), which 
can be calculated using Eqs. (lT5j) -(fT7 |) : 



pqod = 1 s - 1 , Pi +P2 ~ 2 + 2s 
2 B B 2 
(1 + Bf + s 2 - (1 + B)(2s + PlP2 ) 
B 3 



ln(l + B), (31) 



where B — 2pip2 —p\~pi and s = y/1 — piyi — P2- For 
the small-error case (at short time t) this gives 

F^ d « 1 - (pi +pl)/24 - 5pipa/12, Pl , 2 « 1. (32) 

If we use the averaging over the Bloch sphere with 
weight proportional to the probability P nn + P rr of the 
measurement result 0, then 



pqcd 



H 2 + \p\ 2 VT^p~iVT^} 2 + M 2 Pn 



■Pnn *rr 



, (33) 



(the numerator and denominator are averaged sepa- 
rately), which gives 



^ qod _ 2 - pi - p 2 + f P1P2 + \/l - Pi y 7 ! - P2 
av ~ 3[l+pip 2 -(pi+P2)/2] 



(34) 



[Note that instead of using Eqs. ((9^- (fTTj) . the 6-point 
averaging trick can be used separately for the numerator 
and denominator of Eq. ([53"]) - see Appendix.] At short 
times this gives F£° d w 1— (pf +P2) /24— 5piP2/12, same 
as for F£° d [see Eq. 
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FIG. 2: Average state preservation fidelities for the two-qubit 
encoding (compared with no encoding), as functions of the 
one-qubit energy relaxation probability p = 1 — e~ t//Tl (same 
for both qubits, pi = P2 — p)- The solid lines show the QED 
fidelities F^ d and F^ A given by Eqs. pi) and (53)1 . (F^ ed 
assumes averaging over the Bloch sphere with uniform weight, 
while for F^ d the weight is proportional to the probability of 
the "no error" measurement result 0.) The dashed line shows 
the QEC fidelity F^ c , which coincides with F^ n , for which 
the measurement result is ignored, Eqs. (|29l) and (|35|) . The 
dotted line shows the one-qubit fidelity F^ without encoding, 
Eq. (|18p . QEC performs worse than no encoding, while QED 
provides a significant improvement for p < 0.3. 



Figure 2 shows the QED fidelity defined in both ways, 
F^ d and F a q v cd , as functions of the one-qubit relaxation 
probability, assuming similar qubits, p\ = P2 = P- For 
p < 0.3 both fidelities are significantly higher than the 
fidelity F^ for an unencoded single qubit [given by Eq. 
([15])] , which itself is higher than the fidelity when 
the ancilla measurement result is ignored [Eq. ([21?]) ]. 

Now let us discuss whether or not the state of the main 
qubit can be made closer to \ipi n ) using the measurement 
result information, as in quantum error correction. If the 
measurement result is 0, then the qubit state is described 
by the first and last lines of Eq. ([2"T]) . It is rather obvi- 
ous that in this case no unitary operation can improve 
further the average fidelity [for QEC we are interested in 
averaging with the weight proportional to probability - 
see Eq. (|24p] . This statement is rigorously proven in the 
Appendix. So, no correction should be applied for mea- 
surement result 0. (Actually, a non-unitary operation 
due to p artial collapse can increase the fidelity in this 
case 29-32], but we consider only unitary operations, as 
it should be in the usual QEC.) When the measurement 
result is 1, the main qubit is in the state |0) with proba- 
bility -Prn/(-Pm + Pm) or in the state |1) with remaining 
probability P nT /(P ln + P nr ) - see Eq. (|2"7|) . In the case 
Pi = P2 this is the fully mixed state, and any unitary 
operation does not change it. Thus a meaningful error 
correction is impossible, and therefore F^° c — F£= n (see 
Fig. 2). 

Actually, if P2 > Pi, then a slight improvement of fi- 
delity is possible by applying the 7r-pulse (exchanging 
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states |0) and |1)) when the measurement result is 1. 
This makes the resulting state closer to |1) than to |0), 
and correspondingly on average closer to IV'in), because 
the probability of measuring 1 increases with \ip- m ) being 
closer to The optimality of this procedure for mea- 
surement result 1 is proven in the Appendix. It is easy 
to calculate the fidelity change due to the 7r-pulse (the 
easiest way is to average over the 6 initial states and to 
work with unnormalized states - see Appendix) . The re- 
sulting optimal QEC fidelity for the 2-qubit encoding of 
Fig. 1 is 



^pqoc 



Vi 1 -Pi)( l -P2) min(pi,p 2 ) 



(35) 



C. iV-qubit encoding 

We now extend our discussion of the protocol of Fig. [U 
including N—l ancilla qubits. The encoded state is then 
a|0 ) +/3\1 N ) . The state evolution can be unraveled into 
2 N scenarios depending on which qubits relax. However, 
there are 2 Ar_1 measurement results, and each of them 
corresponds to two scenarios. If the main qubit does 
not relax, then the measurement result directly shows 
which ancilla qubits relax (i.e. result 1 indicates the re- 
laxation event), while if the main qubit relaxes, then the 
relaxation scenario is shown by the complement of the 
measurement result (i.e. result indicates relaxation). 

The measurement result (all zeros) indicates that the 
main qubit is either in the state 



i 



P:i 



(36) 



where P nonc = M 2 + |/3| 2 IIj=i (1 — Pj) 1S the probability 
that no qubits relax, or in the state |0) if all qubits relax, 
with the corresponding probability P a n = rX/=iP.r 
For any other measurement result the main qubit is cither 
in state |0) or |1), with easily calculable probabilities of 
the scenarios. For simplicity we assume pj — p below. 

As in the previous subsection, we consider three possi- 
ble ways to proceed: ignore the measurement result, se- 
lect only cases with measurement result (QED), or try 
to improve the fidelity when an error is detected (QEC). 
If the measurement result is ignored, the average fidelity 
(calculated in a similar way as above) is 
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it obviously decreases with increasing number of ancilla 
qubits. 

In quantum error detection we select only cases with 
measurement result 0. Then the state fidelity is 
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and averaging it over the Bloch sphere with uniform 
weight we obtain 
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ln(I + P), (39) 



where B = —l + (l—p) N +p N and S = 2(1 -p) N/2 +p N . 
For N = 2 this equation corresponds to Eq. (f3"Tj) . The 
small-error approximation for N > 3 is 



F2f « 1 - 7VV/24, p«l. 



(40) 



It is interesting to note that this approximation does not 
work for N = 2, for which F^° d w 1 — _p 2 /2, as follows 
from Eq. (|3"2"1) . The reason is that P a ii scales as p N , and 
therefore for N > 3 Eq. (|4T)|) does not have a quadratic 
contribution from the scenario when all qubits relax; the 
infidelity comes only from the difference between |f/ , none) 
and IV'in)- In contrast, for N — 2 the fidelity P a ^ d is fur- 
ther decreased by p 2 /3 due to relaxation of both qubits. 

From Eqs. (gD} and ([32]) we see that the best QED 
fidelity in the small-error case (p <C 1) is achieved by 
the 3-qubit encoding, N = 3; then 1 - F^ d « (3/8) p 2 . 
However, this is only the factor 4/3 better (smaller) than 
for N — 2. Therefore, from the experimental point of 
view the 2-qubit encoding (which is easier to realize than 
the 3-qubit encoding) seems to be most natural. 

If we use the averaging of the QED state fidelity (|3"5)) 
with weights proportional to the probability P non c + Pali 
of the measurement result 0, then we essentially need 
to average the numerator and denominator of Eq. 
separately, thus obtaining 
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(41) 



In the small-error case (p <C 1) for N > 3 this gives 
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1 - N 2 p 2 /24, same as Eq. (gDJ) for F\ 
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Figure 3(a) shows the QED fidelities P a q v od and F^ d 
for N = 2, 3, and 4. We see that the difference between 
F^ d and F^ d becomes larger with increasing N, but 
the difference is small at small p. Note that the QED 
fidelity for the 2-qubit encoding becomes better than for 
the 3-qubit encoding for p > 0.3. 

Now let us discuss the possibility of QEC proto- 
cols, which use unitary correcting operations depend- 
ing on the measurement result. If the result is 0, then 
the unnormalized density matrix of the main qubit is 
Pnono|'0nono)('0nonc| + P a ii|0)(0|. As proven in the Ap- 
pendix, no unitary operation can increase the fidelity 
in this case (in contrast to non-unitary partial-collapse 
operations [29h32| |). For all other measurement results, 
the main qubit is in the incoherent mixture of the 
states |0) and |1); the unnormalized density matrix is 
P),m|0)(0| + Pi, m |l)(l|) where the corresponding proba- 
bilities are Pi, m = |/3| 2 (1 - p) ]\ i=2 N /(m,-) and P , m = 



7 




ft, 

'6- - 



1.0 
0.9 
0.8 
0.7 
0.6 



0.5 
0.0 





N=2 


(b) 














■ /7l eC - 
* av 






F ign -■ 
1 av 


^^^^^^ 




r av 







0.2 



0.4 0.6 
P 



0.8 



1.0 



FIG. 3: (a) The QED fidelities F^ d (solid lines) and F^ d 
(dashed lines) for the encoding using N = 2, 3, and 4 phys- 
ical qubits, as functions of the single-qubit energy relaxation 
probability p. The dotted line shows the fidelity for an 
unencoded qubit. (b) The optimal QEC fidelity F^ c (solid 
lines) and the fidelity F^ n when the measurement result is 
ignored (dashed lines) for N = 2, 3, and 4. 



|/3| 2 ?n i= 2,jv /(1-mO, where /(l) = p, /(0) = 1-p, and 
mi is the measurement result for the ith ancilla qubit. 
As shown in the Appendix, the maximum fidelity is then 
achieved by applying the 7r-pulse (exchanging |0) and |1)) 
if Pi >m < -Po,m and doing nothing if Pi >m > Po,m- Calcu- 
lating the corresponding qubit state fidelity (compared 
with the initial state), summing over the 2 Ar_1 measure- 
ment results, and averaging over the Bloch sphere, we 
obtain the QEC fidelity 
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maxtp-p^, (1-p)- (1-p)*]. (42) 



The QEC fidelity as well as the fidelity F^f 1 for ignor- 
ing the measurement result are shown in Fig. 3(b) for 
N = 2, 3, and 4. The curves for F°£ c and F^f coincide 
at p < 1/2, because in this case the optimal correction is 
no correction. We see that for any N the QEC fidelity 
is smaller than the no-encoding fidelity F^, so the error 
correction by a repetitive code does not protect against 
energy relaxation. 



D. Discussion 

Our results show that the repetitive codes do not work 
for QEC protection against energy relaxation. This is 
because energy relaxation is very different from a bit 
flip (or phase flip or bit-phase flip), for which repetitive 
codes work well. In the language of quantum stabilizer 
codes [2| the event of energy relaxation corresponds to 
the "sum" of two errors: bit flip (X-operation) and bit- 
phase flip (F-operation) - see the Kraus operator A r in 
Eq. ([5]). So, a stabilizer code should be able to protect 
against both of these errors to protect against energy re- 
laxation events. [Actually, a weaker error due to the "no 
relaxation" Kraus operator A r in Eq. (JS|) also requires 
protection against phase flip (Z-operation) errors.] For 
exam ple, the standard 5-qubit and 7-qubit QEC codes 
[E 13. l24j protect against all 3 types of errors (X, Y, Z), 
and therefore can protect against energy relaxation. 

Using the approach of the stabilizer codes and the 
quantum Hamming bound let us calculate the min- 
imum number of qubits N to protect against X and Y 
errors. The Hilbert space of the dimension 2 N can be 
divided into 2 N ~ 1 orthogonal two-dimensional subspaces 
( "copies" of the qubit space) ; these subspaces should be 
able to distinguish the cases with various errors and no 
error. Since the number of possible errors is 2N, we have 
an inequality 2 N ~ 1 > 1 + 2N. From this inequality we 
find N > 5 (the same minimum as for all 3 types of quan- 
tum errors). Notice, however, that an approximate QEC 
for energy relaxation is possible for N = 4 [33[ (see also 
34(|). This code breaks the above limitation because the 
relaxation event is treated as one error, not as a "sum" 
of X and Y (the drawback though is a slightly proba- 
bilistic operation). In any case, the QEC codes protect- 
ing against energy relaxation are much more complicated 
than the repetitive codes. 

Even though the repetitive codes are not good for QEC 
protection against energy relaxation, we have shown that 
they can be well used for QED. Moreover, only 2-qubit 
encoding is sufficient for that. An interesting question 
is whether or not it is beneficial to do many cycles of 
QED, correspondingly decreasing the time of each cy- 
cle and therefore the error probability p in each cycle 
(such division into shorter cycles is beneficial for QEC 
P, HJ). The simple answer is that such division into 
shorter cycles does not help much for the protocol of 
Fig. 1. The reason is that even when no relaxation 
events happen, the qubit state changes - see Eq. (|36p . 
because the absence of relaxation preferentially indicates 
state |0) and plays the same role as the partial collapse 
[30j . Rewriting Eq. (|3^|) in a non- normalized way as 
l^nonc) = a\0) +/3|l)nJ = iexp(-*/2Tij), we see that 
division into several QED cycles does not change the fi- 
nal wavefunction \ip none ) as long as the total time t is the 
same. Therefore, since for N > 3 this evolution is the 
main reason for imperfect QED fidelity at p < 1 (see 
discussion in the previous subsection), there is not much 
benefit of using the QED cycles. Nevertheless, some im- 
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provement of the QED fidelity will be due to a decrease 
of the probability P a n that all qubits relax. Since this 
probability scales as p N « (t/Ti) N , the division into M 
cycles is expected to decrease the corresponding contri- 
bution to the procedure infidelity by the factor M N . 
This improvement will be most significant for N = 2: it 
will essentially change the approximation F^ d ss 1— p 2 /2 
given by Eq. ([32} into i™ cd «l-p 2 /6 given by Eq. CP 
for N = 2. 

A more important improvement of the QED fidelity 
can be achieved if the no-relaxation evolution of 
iV'nono) is compensated. One way is to apply a partial 
measurement [29|-|3l| to the main qubit after the proce- 
dure, essentially eliminating the evolution (|36[) for the 
price of a further decrease of the selection probability 
(probability of success). Another, easier way is to apply 
7r-pulses, exchanging states |0) and |1), between (after) 
the QED cycles (these 7r-rotations can be around any 
axis in the equatorial plane of the Bloch sphere) . Then 
for an even number of equal-duration QED cycles, the 
no-relaxation evolution (|36f will be com pen sated exactly 
(as in the uncollapsing procedure (29l . [3lll32l .l35j). and the 
QED infidelity 1 — F qod will be only due to the contribu- 
tion from P a ii- (The use of 7r-pulses resembles dynamical 
decoupling of the "bang- bang" type [36j]; however, the re- 
semblance is accidental, since dynamical decoupling can- 
not protect against the energy relaxation |37|.) 

For an estimate of the corresponding QED fidelity, let 
us consider the procedure with total duration t < T\ , di- 
vided into M cycles of duration t/M each {M is even). 
In each cycle p « t/MTi <C I, and if we assume 
Nt/MTx < 1, then |V„o„c) « |V*n) in Eq. The 
probability that the iV-qubit relaxation (which remains 
undetected) happened in the first cycle is P a n = \f3\ 2 p N , 
and similar probability for the second cycle (after tt- 
pulse) is lapp^. Therefore, in a selected QED realiza- 
tion (with all measurement results 0) the probability to 
have an undetected relaxation event is (M/2)p N , inde- 
pendent of the initial state (we assume this probability 
to be small, then we can neglect the double-events). If 
such an undetected relaxation event happens, then the 

average fidelity is P av = |a| 2 |/3| 2 /|aF = V 3 - The Q ED 
fidelity then can be calculated as I — (1 — 1/3) x (M/2)p N , 
which gives 

P?v d « P> d « 1 - M(t/MTi) N /3. (43) 



(If the above assumption Nt/MTi <C 1 is violated, then 
the factor 1 /3 changes, but the scaling remains the same.) 
We see that for this procedure the division into a larger 
number of cycles M is beneficial, as well as using more 
qubits (N) for the encoding. Note that our QED proce- 
dure does not prevent the relaxation events from happen- 
ing, so the average probability of observing the "no error" 
result in all M cycles is approximately exp(—tN/2Ti). 
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FIG. 4: Two-qubit experimental protocol for realizing quan- 
tum error detection/correction. Notations E\ and Ei rep- 
resent the four detectable error rotations: Pf (20), Rl(28), 
R Y (26), and i?f(26>). Notations Y/2 and -Y/2 represent 
R Y {n/2) and R Y (-ir/2), respectively. 



III. TWO-QUBIT ERROR DETECTION AND 
CORRECTION FOR PHASE QUBITS 

In this section we propose and analyze the operation of 
two-qubit error detection/correction protocols designed 
for experimental implementation with the current tech- 
nology of superconducting phase qubits. We will discuss 
several similar protocols (including the QED protocol for 
energy relaxation) ; for all of them the goal is to preserve 
an arbitrary initial state \ipi n ) = ct\0) + (3\1) of a qubit. 

The first procedure (which we will mostly consider) 
is shown in Fig. [4j it is designed to preserve the state 
|?/;i n ) of the upper (main) qubit. Encoding is per- 
formed by preparing the lower (ancilla) qubit in the state 
(|0) + |l))/-\/2 by starting with the ground state |0) and 
using the Y- rotation over the angle ir/2 (denoted as Y/2), 
and then applying the controlled- Z (CZ) gate between 
the two qubits. [Note that the CZ gate is the natural 
entangling operation for the phase qubits [19j, 123J, |38j.] 
This produces the entangled two-qubit wavefunction 

[«|0> ® (|0) + ]X» +^|1> ® (|0> - |1»]/V2, (44) 

where the leftmost entry represents the main qubit. Af- 
ter encoding, the decoherence process is simulated by 
applying a unitary rotation to one of the qubits. For 
this encoding we consider a set of four possible rotations: 
Pf {29), RY(29), R Y {29), and Pf (20), where the sub- 
script indicates the qubit number (1 for the main qubit), 
the superscript is the rotation axis on the Bloch sphere, 
and the argument 29 is the rotation angle on the Bloch 
sphere (the corresponding rotation angle in the wavefunc- 
tion language is 9). 

After the error rotation has been applied, the resultant 
state is decoded by inverting the encoding operation, and 
the ancilla qubit is measured in the computational basis. 
In the absence of error rotation, the state after decod- 
ing is (a|0) + /?|1)) (g) |0), so that the initial state of the 
main qubit is restored and the measurement result for 
the ancilla qubit is 0. The error rotation disturbs the 
final state, which probabilistically changes the measure- 
ment result to 1 and also changes the final state of the 
main qubit. 
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A. Analysis of the ideal case 

Let us start with analyzing the effect of the error ro- 
tation R x (29) (X-rotation of the main qubit). It trans- 
forms the encoded state (|4"H) into the state 

[a(cos0|O) -isin0|l)) <g> (|0> + |1)) 
+ i S(-ism6»|0> +cos(9|l>) <E> (|0) - |1))]/V2, (45) 

which after decoding (but before measurement) becomes 

cos0|Vin)|O} + isin0X|V> in }|l}, (46) 

where X is the Pauli- A-matrix transformation . It is 
clear that we obtain ancilla measurement result with 
probability cos 2 9, and then the state of the main qubit 
is restored to l^in), or obtain result 1 with probability 
sin 2 0, which leaves the main qubit in the state X|^j n ). 
[In this section we use the standard quantum comput- 
ing notations [2|, in which the Pauli matrices act on 
column vectors with the upper element corresponding 
to the state |0). Note that for one-qubit wavefunctions 
R x (n) = -iX] 

In quantum error detection we select only result 0, and 
this gives the perfect state preservation fidelity, F s ^ cd = 1, 
for any initial state. We can also use the approach 
of quantum error correction and apply the X gate [i.e. 
R x (it)] to the main qubit when the error result 1 is mea- 
sured. This produces the initial state IV'in) for both mea- 
surement results with perfect fidelity, F^ c = 1. There- 
fore, the QED and QEC fidelities averaged over the Bloch 
sphere are also perfect, 

Fr 1 = Fr = i. (47) 

Notice, however, that for QEC we had to know that an 
error is due to the X-rotation applied to the first qubit. 
This is different from "real" error correction, in which 
we do not know the type of error, but is acceptable for a 
demonstration experiment. 

Let us also calculate the storage fidelity if the mea- 
surement result is ignored (or, equivalently, the an- 
cilla qubit is not measured). From Eq. (|46|) we obtain 
the state fidelity for the main qubit F*f n = cos 2 9 + 
sin 2 9 (ipin\X\ipi n ) 2 , which after averaging over the Bloch 
sphere becomes 

F^ n = cos 2 9 + (sin 2 9) /3. (48) 

Note that if the rotation R x (29) is applied to a qubit 
without encoding, then the average fidelity is still given 
by Eq. (l48l) . so the encoding with ignored measurement 
result (or no measurement) does not affect the average 
preservation fidelity (moreover, it does not affect the 
state fidelity for any initial state). 

Now let us analyze in a similar way the case when the 
error is introduced by the Y-rotation of the main qubit, 
R\ (29). Then the two-qubit state before the measure- 
ment is 

cos0|^)|O)+*sm0Y|V>in>|l), (49) 



so that the measurement result still restores the initial 
state IV^n) of the main qubit, while for the measurement 
result 1 the state of the main qubit is Y\ipi n ), thus re- 
quiring the Y-gate correction [i.e. R Y (ir) — —iY]. The 
QED and QEC fidelities are still perfect, Eq. (J47J), while 
the fidelity with ignored result is still given by Eq. (|48l) . 
Note that the correcting Y-gate is different from the cor- 
recting X-gate in the previous case, so we need to know 
the type of the error to apply the proper correction (in 
a demonstration experiment the error rotation is applied 
intentionally, so its type is obviously known). 

Now let us consider the error due to the Y-rotation 
of the ancilla qubit, R^ (29). Then the state before the 
measurement is 

cos0|V>m)|O} + sin0Z|V>in)|l), (50) 

and therefore in the case of measurement result 1 the Z- 
gate correction is needed to restore IV'in), while for the 
measurement result no correction is needed. Equations 
(1471) and (gHJ) are still valid. 

Finally, for the Z-rotation of the ancilla qubit, i?f (20), 
the state before the measurement is 

|V>in)(cos0|O)+isin0|l)). (51) 

The final state of the main qubit is insensitive to this 
rotation, and therefore no correction is needed for both 
measurement results. In this case Eq. (|4"T|) is still valid, 
while Eq. (g5) is replaced by i™ n = 1. 

We have discussed the effect of four error rotations: 
R x (29), Rj(29), R%(26), and i?f (29). The two re- 
maining rotations, i?f (20) and R x (29) 7 gradually change 
the final state of the main qubit (both produce its Z- 
rotations) but always produce final state |0) of the ancilla 
qubit. Therefore, these errors are undetectable and are 
excluded from our set of error rotations. 

As discussed above, for QEC we need to know which 
one out of four error types has been applied. In contrast, 
for QED we do not need to know the error type; for all of 
them the measurement result indicates the perfect state 
of the main qubit. Moreover, for QED these types of error 
rotations can be applied simultaneously, as long as the 
rotation angles are relatively small, to make negligible the 
second-order terms (which in the QEC/QED language 
correspond to double-errors). 

It is most natural to view the analyzed procedure as a 
QED protocol. However, we would like to emphasize that 
its interpretation as a QEC protocol is also possible: the 
proper correction is possible when we know which error 
process is applied. (In the existing QEC experiments the 
types of allowed errors are almost always limited; in our 
protocol the number of allowed types is further reduced 
to one out of four.) Most importantly, our simple two- 
qubit protocol demonstrates the main "miracle" of QEC, 
that continuous quantum errors can be transformed into 
discrete errors, and then corrected. 
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B. Realization using phase qubits 

So far we considered the ideal case when there is no 
physical decoherence of qubits, and the loss of fidelity is 
only due to intentional rotations of the qubit states. In 
this subsection we discuss a more realistic experimental 
situation, with added decoherence during the protocol. 
We will have in mind the present-day technology of su- 
perconducting phase qubits 0, [H, |38| . 

Note that the phase qubit techn olog y provides a high- 
fidelity measurement (about 95% pjj], so we consider it 
perfect in the simulations); however, it takes a significant 
time to read out the measurement result (longer than the 
qubit decoherence time) . While this is not a problem for 
the QED, the QEC at present cannot be done in real 
time. Nevertheless, there is a simple way to go around 
this difficulty in an experiment. The resulting state of 
the main qubit is measured by using the quantum state 
tomography (QST), so the experiment is necessarily re- 
peated many times. It is easy to separate the QST data 
for ancilla measurement results and 1. In this way 
two different density matrices of the main qubit are ob- 
tained for ancilla measurement and 1. For the result 1 
it is then easy to calculate the density matrix after the 
correcting operation (if it were applied in real time) . Fi- 
nally adding the two density matrices (with weights equal 
to the probabilities of ancilla measurement results), the 
qubit density matrix for the QEC procedure is obtained. 

In phase qubits [ltl [H, |38[ the main sources of deco- 
herence are single-qubit energy relaxation (with T± on the 
order of 0.5 /lis) and pure dephasing (with a comparable 
or a little shorter dephasing time T v ). The decoherence is 
somewhat reduced in the RezQu architecture [H, [2(| , in 
which the quantum information is often transferred be- 
tween the phase qubits and resonators (resonators have 
much longer T\ and practically no pure dephasing). We 
have simulated the procedure of Fig. 4 in a simplified 
way, which does not explicitly reproduce the RezQu im- 
plementation of the protocol, but still uses a reasonable 
account of realistic decoherence. 

For simplicity for each qubit we assume T\ = Ti (so 
that the pure dephasing time is T v = 2Ti). We assume 
that single-qubit rotations [including R Y (±ir/2) of an- 
cilla qubit, preparation of the main qubit, and error ro- 
tations] take 10 ns each, CZ gates take 40 ns each, and 
there are 5 ns spacings between the operations. Then the 
whole protocol of Fig. 4 (ending before measurement of 
ancilla qubit and tomography of the main qubit) takes 
135 ns. We calculate the evolution of the two-qubit den- 
sity matrix by breaking the procedure into small time 
steps and applying energy relaxation and pure dephasing 
to each qubit (for simplicity the CZ gate is simulated as 
a gradual accumulation of the phase, as would be for the 
dispersive gate). We start with 6 initial states of the main 
qubit [|0), |1), (|0)±|l))/ v / 2, (|0)±z|l))/v / 2, labeled by 
index j = 1,...6 below], and from the final two-qubit 
density matrices we calculate reduced unnormalized one- 
qubit density matrices poj and pij, corresponding to 
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FIG. 5: Numerical results for the average QED fidelity F^ d 
(solid lines), the QEC fidelity f^ v cc (dottes lines), and the 
fidelity Fl^ n with ignored ancilla measurement results (dashed 
lines), as functions of the angle 28 of intentional X-rotation of 
the main qubit, (28). The simulated protocol of Fig. 4 has 
duration of 135 ns. We assume the qubits with Ti = Ti — 300 
ns, 500 ns, and 700 ns. 



ancilla measurement results and 1 (the probabilities 
of these results are then Trpoj and Txpij). The aver- 
aged preservation fidelity with ignored measurement re- 
sults is then (see Appendix) = (1/6) Y\ Tr[(p 0ij + 

/Oi,j)p™], where p™ — \ipj){ipj\ is the unchanged initial 
state. The averaged (weighted) QED fidelity is then 
^v° d = V, I '•:/<;.. ,/' H ' : ' : i V/lr/A,,. and the QEC fi- 
delity is F£ c = (1/6) £\ Tr[(p 0j + PlJ)pf], where the 
corrected density matrix p1°" is obtained from pij by 
applying the ideal correcting operations (X, Y, Z, I) dis- 
cussed in the previous subsection. 

Figure 5 shows the average fidelities F^ d (solid lines), 
F^ c (dotted lines), and F^f 1 (dashed lines), as functions 
of the rotation angle 28 (in units of 7r) for the intentional 
X-rotation of the main qubit, Rf(26). The three sets 
of lines are for three values of T% = T2: 300 ns, 500 ns, 
and 700 ns. Note that we present the average fidelities 
F avi but they can be easily converted into the process 
matrix fidelities F x via Eq. @ . Also note that the range 
from 1/3 to 1 for F av (used for the vertical axis in Fig. 
5) corresponds to the range from to 1 for F x . 

From Fig. 5 we see that even for T\ = T% = 300 ns the 
QED fidelity is significantly higher than the fidelity with 
ignored measurement result (recall that the procedure 
duration is 135 ns). The difference between F^° d and 
i^v" becomes larger for longer decoherence time (500 ns 
and 700 ns). The QEC fidelity is below the QED fidelity 
(and even below i^v") f° r small 9, but becomes above 
^ v ed and Flf 1 at large 6. 

It is interesting to notice that F^| n at 20 s» it is much 
closer to the ideal value 1/3 [see Eq. ([IS])] than to the 
ideal value 1 at 29 « 0. This property can be under- 
stood using the equivalent language of the process fidelity 
F x = Tr(xdcsircdX) _ see Eq- ©■ Since the desired opera- 
tion is the absence of evolution, F x = xn in the standard 
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FIG. 6: Same as in Fig. 5, but for four types of intentional 
qubit state rotations: R?(26), R\(26), Rt(26), and i?f (29). 
The qubits with Ti = T2 = 500 ns are assumed. Results for 
rotations R* (29) and RX (29) practically coincide. 
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FIG. 7: Modified two-qubit QED/QEC algorithms. The pro- 
tocol shown in (a) detects/corrects errors due to rotations R{ , 
R\ 1 fl^T, and i?f ; it can be used to protect from natural pure 
dephasing of the qubits. The protocol in (b) is designed for 
error rotations R\ , R\ , R2 , and R\ . Therefore, it can be 
used as a QED procedure for errors due to energy relaxation 
of both qubits (stored in resonators) . 



C. Related protocols 



notations for the one-qubit 4x4 matrix \ @, E3, HH, HH ; 
note that \n + Xxx + Xyy + Xzz = 1- Ideally \n = 1 
for 29 = and xxx = 1 for 29 = it. Since decoherence 
spreads these ideal unity values to the three other diag- 
onal elements of %, we would expect that F^ n at 28 = tt 
should be (very crudely) three times less than 1 — F^ gn 
at 29 — 0. This roughly corresponds to what we see in 
Fig. 5. 

Figure 5 shows the results only for the X-rotation of 
the main qubit, R K (26). The results for all four consid- 
ered rotations, i?f (26), R\(29), Rj(26) and i?f (26), are 
shown in Fig. 6 for T\ — T2 — 500 ns. The results for X 
and Y-rotation of the main qubit are practically indistin- 
guishable from each other. The QED and QEC fidelities 
for Y-rotation of the ancilla qubit are very close to the 
corresponding fidelities for the rotation of the main qubit. 
For i?-rotation of the ancilla qubit the operation with ig- 
nored measurement coincides with the QEC operation 
(because no correction is applied for measurement result 
1), and the QED fidelity 7^ v ed is higher than i^| n = F^ v oc 
only at 26 < tt/2, and only by a small amount. Obvi- 
ously, the rotation i?f (2#) is not good for demonstrating 
an advantage of this encoding, in contrast to other rota- 
tions. 

Overall, from Figs. 5 and 6 we see that the current tech- 
nology of phase qubits is good enough for demonstrat- 
ing the operation of the considered two-qubit QED/QEC 
protocol. In an experiment, the larger value of the QED 
fidelity in comparison with the case of ignored measure- 
ment result is the demonstration that the QED procedure 
is beneficial. Similarly, the QEC operation can also be 
demonstrated (though with the caveat discussed in the 
previous subsection). 



The protocol of Fig. 4 can be easily modified to change 
the set of four detectable/correctable error operations. 
For example, if we desire protection from Y and Z rota- 
tions of both qubits [i.e. Rj(26), Rf(26), R%(20), and 
R 2 (29)}, we can add ±7r/2 Y-rotations of the main qubit 
before and after the error rotations - see Fig. [7[a). Such 
encoding also protects from natural pure dephasing of 
both qubits. 

For protection from X and Y rotations of both qubits 
[i.e. Rf(26), R\{29), R 2 K (29), and R%(26)], we can add 
T-7r/2 ^-rotations of the ancilla qubit before and after 
the errors - see Fig. EJb). Such encoding can be used in 
the QED mode for the energy relaxation of both qubits. 
(This procedure essentially realizes the idea of Fig. 1 for 
two qubits; the only difference is the encoding cv|00) — 
/3 1 1 1 ) instead of encoding a|00) + /3|11) considered in Sec. 
II B.) 

In the RezQu architecture based on phase qubits 
HjJ [26|, the protocol of Fig. [7£b) can be efficiently used 
for storing the information in the resonators. We have 
simulated the operation of this protocol, assuming that 
the encoding/decoding is done with the phase qubits hav- 
ing relaxation times T\ = T 2 , while in between encoding 
and decoding the information is moved to resonators for 
a relatively long storage. The procedure (without the 
storage time) is slightly longer than what was consid- 
ered in the previous subsection: 155 ns instead of 135 
ns (we still do not simulate explicitly the move opera- 
tions between the qubits and resonators). Solid lines in 
Fig. 8 show the corresponding QED fidelities F^ d as 
functions of the single-qubit energy relaxation probabil- 
ity p — 1 — exp ( — £ storage /T 1 rcsonator ) during the storage 
(in experiment [2l] x , 1 resonator > T x , though our results do 
not need this assumption). The QEC operation is impos- 
sible in this case (for real energy relaxation in resonators) ; 
however, as seen from Fig. 8, the QED operation can be 
reliably demonstrated experimentally. 
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FIG. 8: Average QED fidelity F^ d (solid lines) for the two- 
qubit protocol of Fig. [7{b), as a function of the single-qubit 
energy relaxation probability p — 1 — exp(— i s toragc/TT osonat ° r ) 
during information storage in resonators. Dashed lines show 
average fidelity F^ n when the measurement result is ignored. 
The encoding/decoding is done with phase qubits having sig- 
nificantly shorter relaxation times Ti = T2 (300 ns, 500 ns, 
and 700 ns); the assumed duration of the procedure (exclud- 
ing storage time) is 155 ns. 



IV. CONCLUSION 

In Sec. II we have analyzed the performance of N- 
qubit repetitive quantum codes in the presence of energy 
relaxation. As expected, these codes are not usable for 
quantum error correction. However, they can be used for 
quantum error detection. The best QED performance for 
weak energy relaxation is provided by the 3-qubit repet- 
itive encoding [see Eq. (|40|) and Fig. 3(a)], while the 2- 
qubit encoding is sufficient and gives only slightly lower 
fidelity [sec Eq. (|3"2"1) ]. We have found that the main con- 
tribution to the QED infidelity for A^ > 3 comes from 
the non-unitary change of the quantum state in the case 
when no relaxation happens. Therefore, the QED infi- 
delity can be strongly decreased if the QED algorithm 
is complemented with partial quantum measurement or, 
alternatively, if the protocol is divided into the even num- 
ber of cycles and complemented with 7r-pulses in between 
(this resembles dynamical decoupling, though only super- 
ficially). In this case the fidelity improves with dividing 
the total duration into a larger number of cycles and us- 
ing more qubits for the encoding [see Eq. (|43|)]. 

Note that the QED fidelity cannot be introduced in 
the usual way [as F x = Tr(xdesiredX)] v i a the process x~ 
matrix. In the analysis we have used the state fidelity, 
averaged over the Bloch sphere, with the usual conver- 
sion into F x via Eq. (HI). For the state fidelity averaged 
with the weight proportional to the selection probability 
(denoted F^° d ) the usual trick of averaging over the 6 ini- 
tial states can be used (see Appendix), while for the state 
fidelity averaged with uniform weig ht (denoted F a q v cd ) we 
had to average over the Bloch sphere explicitly, using 
Eqs. (|5))- (fT7|) . In the analysis we used unraveling of the 
decoherence evolution into the "error scenarios" , which 



is in general similar to the standard approach used in the 
quantum error correction, but is different in the way that 
each scenario describes a non-unitary process. 

In Sec. Ill we have considered simple two-qubit pro- 
tocols of quantum error detecton/correction, suitable 
for present-day experiments with superconducting phase 
qubits 23]. In the protocol of Fig. 4 the errors are simu- 
lated by intentional unitary rotations of the qubit states 
(two types of rotations for each qubit). In this case not 
only the QED, but also the QEC operation is possible 
if we know the applied type of error rotation. Most im- 
portantly, this experiment would demonstrate the QEC 
"miracle" of converting continuous quantum errors into 
discrete errors, which are then correctable. The numer- 
ical simulations (Figs. 5 and 6) with account of deco- 
herence during the protocol show that the experimental 
QED and QEC fidelities are expected to be significantly 
higher than the fidelity with ignored result of the ancilla 
qubit measurement. Therefore, the QED and QEC ben- 
efits can be demonstrated experimentally. 

A slightly different protocol, shown in Fig. 7(b), can 
be used as a QED procedure for errors due to natu- 
ral energy relaxation of qubits stored in resonators of a 
RezQu-architecture device [H,|2^] based on phase qubits. 
The numerical simulations (Fig. 8) show that such ex- 
periment can also be realized with the present-day tech- 
nology, demonstrating the benefits of encoding a logical 
qubit in several (two in this case) physical qubits. While 
the measurement-free QEC experiment has been recently 
realized with superconducting transmon qubits [15[, the 
experiments proposed and analyzed in this paper would 
be the first measurement-based QED/QEC protocols re- 
alized with superconducting qubits. 

The authors thank John Martinis, Andrew Cleland, 
and Matteo Mariantoni for fruitful discussions. The work 
was supported by the IARPA under ARO Grant No. 
W911NF-10-1-0334. 



Appendix A 

In this appendix we prove that no additional unitary 
operation can improve the fidelity of the protocol dis- 
cussed in Sec. II (for 2-qubit or iV-qubit encoding) in 
the case of "no error" measurement result (0 or 0). We 
also prove that for a measurement result which indicates 
an error, the optimal correction is either identity or the 



7r-pulse, exchanging 0) and |1). Along the way we also 
discuss the trick 27], 28| of using only 6 initial states for 
averaging the state fidelity over the Bloch sphere. 

Let us first consider an arbitrary (not necessarily trace- 
preserving) linear one-qubit quantum operation, which 
transforms initial states |1), |0), |±) = (|0) ± |1))/V2, 
I ± i) = (|0) ± i\l))/V2 into the density matrices po, 
Pi: P±i P±i- The center of the Bloch sphere (|0)(0| + 
|l)(l|)/2 = 1/2 is transformed into p c . Because of the 
linearity, only four linearly independent initial states are 
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sufficient to define the operation. So, for an initial state 
with the Bloch sphere coordinates {x, y, z}, 

Pin = (I + xX + yY + zZ)/2, (Al) 

where {X, Y, Z} are the Pauli matrices (x = 1 corre- 
sponds to |+), y = 1 corresponds to | + i), z = 1 corre- 
sponds to |0)), the final state is 

Pfin = Pc + x(p + - p c ) + y(p+i - p c ) + z(po - Pc)- (A2) 

To compare this operation with a unitary U , we calcu- 
late the state fidelity Tr(jOfi n pg n ) (the superscript U in a 
notation means that it relates to the unitary U), 

F st = Tr{[p c + x{p + - p c ) + y(p +l - p c ) + z(p a - p c )] 
x [p u c + x{p u + p u c ) + y{p u +l p u c ) + z{pV - p u c )]}. (A3) 

Note that p v c = 1/2, since a unitary operation does not 
change the Bloch sphere center. In averaging F st over the 
Bloch sphere we average over the coordinates {x, y, z} 
and use the obvious relations x = y = 1 = xy — ~xz — 
yz = 0, x 2 = y 2 = z 2 — 1/3, this obtaining 

F = \^ P c + ^Tr[(p+ - Pc ){p u + - p u c ) 

+ ( P+l - p c ){p u +i -p u c ) + { PQ - p c )( p u -p u c )]. (A4) 

Note that in general we deal here with non-normalized 
density matrices, in contrast to the formalism used in 
Sec. II. Therefore, compared with notations of Sec. II, 
F = F av = F av only for a trace-preserving operation, 
while for a non-trace-preserving operation F = F av P, 
where P is the average probability of selection [see Eq. 
(|24p]. and there is no direct relation between F and F av . 

Using Eq. (|A4[) it is easy to see why averaging over 
the Bloch sphere is equivalent to averaging over only 6 
initial states: |0), |1), |±), and | ± i). The state fidelity 
F + for the initial state |+) is given by Eq. (|A3[) with 
x = 1 and y = z = 0. The state fidelity F_ for the initial 
state |— ) is similar, but x = — 1. It is easy to obtain the 
sum, F+ + F- = Tr p c + 2Tr[(p+ - p c )(p l ( - p")], which 
is similar to the terms in Eq. (|A4[) . Similarly finding the 
sums F +i + F^i and Fq + F\, we obtain [23, [28( 

F = (F + F t + F+ + F_ + F +i + F_i)/6. (A5) 

Note that this relation remains valid for non-trace- 
preserving operations, when we are working with a lin- 
ear operation and non-normalized states. The same six- 
point-averaging relation is valid for the average proba- 
bility of selection P, because P — Trpg n and therefore 
P = Trp c (even two-point averaging is sufficient for P, 
when we choose two opposite points on the Bloch sphere). 
Therefore the six-point-averaging trick is useful for find- 
ing F av = F/P. 

Now let us discuss why an additional unitary cannot 
improve the QEC protocols of Sec. II when the "no error" 
measurement result (or 0) is obtained. The final state 
in this case is an incoherent mixture of the results of two 
linear operations: 

a\0)+j3\l) ->a|0)+fy3|l), a|0) +j0|l) -> k0\Q), (A6) 



where the real positive constants fc and k should obvi- 
ously satisfy inequality k 2 + k 2 < 1. For this operation 
it is easy to find explicitly 

Pc = (1 + k 2 + k 2 )I/4 +(l + ~k 2 - k 2 )Z/4, (A7) 
p+ = Pc + kX/2, p +t = p c + kY/2, po = (Z + I)j2. (A8) 

Then using Eq. (|A4[) we obtain 

F = ~(1 + k 2 + k 2 ) + ~Tr[kX( p v - P u c ) 

+kY(pV i - P V) + l - k2 2 + k2 Z{p u - p u c )\- (A9) 

(Note that comparing the operation with U we assume 
the correction operation U'.) Optimizing each term un- 
der the trace over the unitary U separately, we see that 
the first term is maximized by unitaries, which trans- 
form |+) — > |+); the maximum for the second term is 
achieved when | + i) — > \ + i), and the maximum for the 
third term is achieved when |0) — > |0) (note that k > 
and 1 — k 2 + k 2 > 0). Since the no-evolution operation 
satisfies all these conditions, it provides the maximum 
fidelity, 

C b „, = / , y_l±*±*±PZ?. (A10) 

Note that the average probability of the process (|A6[) 
is P = (1 + k 2 + k 2 )/2, so F av = (2/3)(l + k + k 2 + 
k )/(l + k 2 + k 2 ). In particular, this is an alternative 
way of deriving Eq. (|34|) by using k = i/l — PiV 1 — P2 
and k = ^pxpz- 

Now let us discuss what is the optimal unitary correc- 
tion operation after obtaining the measurement result 1 
in 2-qubit encoding or any result except in iV-qubit en- 
coding. Then the resulting state is an incoherent mixture 
of two linear operations: 

a|0)+/3|l) ->fcj0|l), a|0) -> fc/3|0). (All) 

Finding explicitly the resulting density matrices 

Po =0, p c = p + = p +i = y|0)<0| + y |1)(1|, (A12) 
we obtain from Eq. (IA4[) 

1,2 ,1,2 u2 _ p 

F = + k -^f- Tr[Z(pV - pV)}. (A13) 

Therefore if A: > k, then the maximum fidelity F max = 
(2k 2 + k 2 )/6 is achieved for any unitary U, which does 
not change |0) (same for the correcting operation W). 
However, if k < k, then the optimal U transforms |0) — > 
|1) (same for W) and F max = (k 2 + 2fc 2 )/6. 
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